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Theoretical  investigations  of  the  electronic  structure  and  radiative 
transition  probabilities  of  atmospheric  spec  ies  and  of  kinetic  reaction  rates 
for  electron-ion  ion-atom,  and  ion-molecule  reactions  are  in  progress. 

These  investigations  are  being  carried  out  using  quantum  mechanical  methods 
and  digital  computer  codes  which  represent  state-of-the-art  numerical 
techniques.  The  electronic  structure  and  static  molecular  properties  are 
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*• being  studied  using  initio  configuration interaction  wavefunctions  which 
have  sufficient  flexibility  to  properly  connect  with  the  ground  and  excited 
asymptotic  atomic  and  ionic  limits.  The  charge  transfer  cross-sectionsare 
computed  using  a  semi-classical  close-coupling  code  which  has  been  developed 
in  this  Center.  Modified  Nutnerov  integration  techniques  are  employed  to 
evaluate  the  matrix  elements  for  the  dissociative-recombination  cross- 
sections.  The  branching  into  the  reaction  product  channels  can  also  be 
calculated  using  a  full  quantum  multichannel  close-coupling  code.^ 

Our  previous  theoretical  study  (AFCRL-TR-75-0509)  of  the  recombination 
kinetics  of  e  +  N0+  had  been  directed  mainly  toward  the  calculation  of  the 
energy  temperature)  dependence  of  the  recombination  cross-section.  The 
pertir.-nt  potential  energy  curves  were  calculated  using  ab  initio  Cl  methods 
and  were  defined  to  near-spectroscopic  accuracy.  An  important  problem  is 
the  determination  of  the  branching  into  NC^S)  and  N(^D) ,  which  has  been 
shown  by  our  studies  to  be  governed  by  states  of  ^Y+  symmetry  of  NO.  A 
calculation  of  the  detailed  branching  from  an  initial  recombination  into  the 
I  state  of  NO  to  atomic  separations  of  N(^S)  +  0(^P)  and  N(^D)  +  0(^P) 
was  the  principal  thrust  of  our  CY79  research  effort. 

Calculations  of  the  branching  ratio  between  the  AZIT  and  1  Y  molecular 
statesof  NO  have  been  carried  out  for  the  dissociative-recombination  of 
e  +  N0+.  The  dominant  recombination  channels  for  this  reaction  for  incident 
energies  less  than  'v,  0.5  eV,  are  the  B“fl,  L^Jl  and  B1  states  of  NO.  How¬ 
ever,  "v  20%  of  the  recombination  occurs  along  the  ^T.  symmetry  path.  Our 
calculated  diabatic  potential  curves  for  ^Y-+  symmetry  indicate  a  maximum 
Hamiltonian  interaction  term,  ,  of  270  cm-1  for  an  internuclear  seperation 
of  2.842  bohrs.  Using  a  full  close-coupling  code,  we  calculate  that  the 

coupling  cross-section  between  the  I  ^Y+  and  A  ^Y+  states  of  NO  never 
-20  2 

exceeds  10  cm  for  collisional  energies  <^2.0  eV.  The  recombination 
along  the  I  channel  thus  follows  a  diabatic  reaction  path  and  we  pre- 

diet  exclusively  that  N[ZD]  +  0  [JP]  are  the  atomic  products  of  dissociative- 
recombination  of  e  +  N0+. 

Preliminary  studies  of  the  0+  +  N2  reaction  were  carried  out  t(£  examine 
likely  reaction  pathways.  The  molecular  correlation  diagram  of  N2O  has 
been  constructed  for  collision  energies  up  to  6.16  eV,  and  including  all 
channels  below  0+[^S]  +  N2  ]•  Spin  correlations  indicate  that  the 

lowest  adiabatic  surface  should  be  of  quartet  symmetry.  Previous  studies, 
using  SCF  procedures,  produced  a  co-linear  reaction  path  (^Y  )  which 
exhibited  a  barrier  height  of  approximately  7  eV.  In  contrast, the  best 
experimental  estimate  of  the  barrier  in  this  system  is  'v  0.2  eV. 

-V  ur  preliminary  studies  indicate  that  the  spin-recoupling  that  occurs 
in  going  from  reactants  to  products  is  not  adequately  represented  with  SCF 
or  MC-SCF  frameworks.  Detailed  calculations  are  required  to  estimate  the 
spin-correlation  errors  and  such  studies  are  proposed  as  a  logical  extension 
of  our  present  research  efforts. ^ 
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The  release  of  certain,  chemical  species  into  the  upper  atmosphere 
results  in  luminous  clouds  that  display  the  resonance  electronic-vibrational- 
rotation  spectrum  of  the  released  species.  Such  spectra  are  seen  in  rocket 
releases  of  chemicals  for  upper  atmospheric  studies  and  upon  reentry  into 
the  atmosphere  of  artificial  satellites  and  missiles.  Of  particular  interest 
in  this  connection  is  the  observed  spectra  of  certain  metallic  oxides  and  air 
diatomic  and  triatomic  species.  From  band  intensity  distribution  of  the 
spectra,  and  knowledge  of  the  f-values  of  electronic  and  vibrational  transi¬ 
tions,  the  local  conditions  of  the  atmosphere  can  be  determined  (Ref.  1). 

Present  theoretical  efforts  which  are  directed  toward  a  more  complete 
and  realistic  analysis  of  the  transport  equations  governing  atmospheric 
relaxation,  including  chemical  effects,  and  the  propagation  of  artificial 
disturbances  require  detailed  information  on  atomic  and  ionic  reaction 
rates  and  on  thermal  opacities  and  LWIR  absorption  in  regions  of  temperature 
and  pressure  where  molecular  effects  are  important  (Refs.  2  and  3).  Although 
various  experimental  techniques  have  been  employed  for  both  atomic  and 
molecular  systems,  theoretical  studes  have  been  largely  confined  to  an 
analysis  of  the  properties  (bound-bound,  bound-free  and  free-free)  of  atomic 
systems  (Refs.  4  and  5).  This  has  been  due  in  large  part  to  the  unavail¬ 
ability  of  reliable  wavef unctions  for  diatomic  molecular  systems,  and 
particularly  for  excited  states  or  states  of  open-shell  structures.  More 
recently,  (Refs.  6-9)  reliable  theoretical  procedures  have  been  prescribed 
for  such  systems  which  have  resulted  in  the  development  of  practical 
computational  programs. 

The  theoretical  analysis  of  atmospheric  reactions  requires  the  knowledge 
of  the  electronic  structure  of  atoms,  ions  and  small  molecular  clusters  of 
nitrogen  and  oxygen  and  the  interaction  of  water  or  other  small  molecules 
with  these  clusters.  Knowledge  of  the  chemistry  of  metal  oxide  species, 
which  might  be  present  in  a  contaminated  atmosphere,  is  also  desired.  In 
addition  the  basic  collisional  processes  involving  electrons,  ions  and 
neutral  particles  must  be  understood  to  evaluate  the  dynamic  effects  of  the 
chemistry  of  the  atmosphere.  One  important  reaction  is  charge  transfer 
in  N+  +  0  ->  N  +  0+.  This  has  been  studied  experimentally  by  Neynaber  (Refs.  10 
and  71)  at  collision  energies  between  0.5  eV  and  25  eV  but  low  energy  (<  0.5 
eV)  data  are  apparently  extremely  difficult  to  measure. 

Because  of  the  difficulty  of  conducting  experiments  to  measure  the 
appropriate  cross  sections  for  many  thermal  energy  atmospheric  processes, 
the  development  of  a  sound  theoretical  method  for  calculating  low-energy 
cross  sections  appears  necessary.  Although  relatively  little  work  of  this 
nature  has  been  done  in  the  past,  enough  theoretical  work  is  available  to 
indicate  that  the  development  of  such  procedures  can  be  made  practical, 
particularly  if  good  wave  functions  and  potential  energy  curves  are  avail¬ 
able  for  the  interacting  species. 
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Earlier  research  (Technical  Report  No.  AEWL-TR-72-1)  was  undertaken  to 
assess  the  reliability  of  theoretically  predicting  diatomic  transition  proba¬ 
bilities.  The  systems  studied  under  this  program  included  the  nitrogen  first 
positive  system,  the  oxygen  Schumann-Runge  system  and  the  nitric  oxide  beta 
system.  In  addition,  a  preliminary  study  of  theoretical  techniques  was 
undertaken  for  predicting  cross  sections  for  collisions  involving  electrons 
and  molecular  ions.  Large  two-body  electron-ion  recombination  coefficients 
have  been  observed  in  recombining  plasmas,  both  in  the  laboratory  and  in 
the  upper  atmosphere.  Bates  was  the  first  to  suggest  that  these  large 
recombination  coefficients  could  be  explained  in  terms  of  the  recombination 
of  electrons  and  molecular  ions  in  a  process  called  dissociative  recombination; 
whereas,  the  coefficients  associated  with  the  recombination  of  electrons  and 
atomic  ions  were  expected  to  be  several  orders  of  magnitude  smaller. 

As  an  extension  of  our  collisional  studies  of  fT1"  +  0,  an  examination  of 
the  branching  ratios  in  the  dissociative-recombination  of  e  +  NO  N[^D,  aS] 

+  0  was  carried  out.  This  study  utilized  the  results  of  our  previous 
research  of  the  electronic  structure  of  N0+  and  NO.  In  addition,  a  preliminary 
examination  of  the  ion-molecule  reation,  0+  +  N2  -*■  N0+  +  N  was  carried  out 
to  access  possible  reaction  pathways.  The  goal  of  this  entire  research 
program  was  to  develop  technical  information  concerning  these  systems  which 
is  relevant  to  DNA  interests  in  upper  atmospheric  reactions. 


CURRENT  STATUS  OF  QUANTUM  MECHANICAL  METHODS  FOR  DIATOMIC  SYSTEMS 


The  application  of  quantum  mechanical  methods  to  the  prediction  of 
electronic  structure  has  yielded  much  detailed  information  about  atomic  and 
molecular  properties  (Ref.  7).  Particularly  in  the  past  few  years,  the  avail¬ 
ability  of  high-speed  computers  with  large  storage  capacities  has  made  it 
possible  to  examine  both  atomic  and  molecular  systems  using  an  ah  initio 
approach,  wherein  no  empirical  parameters  are  employed  (Ref.  11).  Ab  initio 
calculations  for  diatomic  molecules  employ  a  Hamiltonian  based  on  the  non- 
relativistic  electrostatic  interaction  of  the  nuclei  and  electrons,  and  a 
wavefunction  formed  bv  antisymmetrizing  a  suitable  many-electron  function  of 
spatial  and  spin  coordinates.  For  most  applications  it  is  also  necessary  that 
the  wavefunction  represent  a  particular  spin  eigenstate  and  that  it  have 
appropriate  geometrical  symmetry.  Nearly  all  the  calculations  performed  to 
date  are  based  on  the  use  of  one-elecrron  orbitals  and  are  of  two  types: 
Hartree-Fock  or  configuration  interaction  (Ref.  8). 

Hartree-Fock  calculations  are  based  on  a  single  assignment  of  electrons 
to  spatial  orbitals,  following  which  the  spatial  orbitals  are  optimized, 
usually  subject  to  certain  restrictions.  Almost  all  Hartree-Fock  calculations 
have  been  subject  to  the  assumption  that  the  diatomic  spatial  orbitals  are 
all  doubly  occupied,  as  nearly  as  possible,  and  are  all  of  definite  geometri¬ 
cal  symmetry.  These  restrictions  define  the  conventional,  or  restricted, 
Hartree-Fock  (RHF)  method  (Refs.  12  and  13).  RHF  calculations  can  be  made 
with  relatively  large  Slater-type  orbital  (STO)  basis  sets  for  diatomic  mole¬ 
cules  with  first  or  second-row  atoms,  and  the  results  are  convergent  in  the 
sense  that  they  are  insensitive  to  basis  enlargement.  The  RHF  model  is  ade¬ 
quate  to  give  a  qualitatively  correct  description  of  the  electron  interaction 
in  many  systems,  and  in  favorable  cases  can  yield  equilibrium  interatomic 
separations  and  force  constants.  However,  the  double-occupancy  restriction 
makes  the  RHF  method  inappropriate  in  a  number  of  circumstances  of  practical 
interest.  In  particular,  it  cannot  provide  potential  curves  for  molecules 
dissociating  into  odd-electron  atoms  (e.g.,  NO  at  large  internuclear  separa¬ 
tion),  or  into  atoms  having  less  electron  pairing  that  the  original  molecule 
[e.g.,  O2  -*■  0(JP)];  it  cannot  handle  excited  states  having  unpaired 
electrons  (e.g.,  the  states  of  O2  responsible  for  the  Schumann-Runge 
bands);  and,  in  general,  it  gives  misleading  results  for  molecules  in  which 
the  extent  of  electron  correlation  changes  with  internuclear  separation. 

Configuration- interaction  (Cl)  methods  have  the  capability  of  avoiding 
the  limitations  of  the  RHF  calculations.  If  configurations  not  restricted  to 
doubly-occupied  orbitals  are  included,  a  Cl  can,  in  principle,  converge  to 
an  exact  wave- function  for  the  customary  Hamiltonian.  However,  many  Cl  cal¬ 
culations  have  in  fact  been  based  on  a  restriction  to  doubly-occupied 
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orbitals  and  therefore  retain  many  of  the  disadvantages  of  the  RHF  method 
(Ref.  8).  The  use  of  general  Cl  formulations  involves  three  considerations, 
all  of  which  have  been  sat isf actor ilv  investigated:  the  choice  of  basis 
orbitals,  the  choice  of  configurations  (sets  of  orbital  assignments),  and 
the  specific  calculations  needed  to  make  wavefunctions  describing  pure  spin 
states  (Ref.  6).  The  first  consideration  is  the  art  associated  with  quantum 
mechanical  electronic  structure  calculations.  Many  methods  (iterative  NSO, 
perturbation  selection,  first  order  Cl,  etc.)  have  been  advocated  for  the 
optimum  choice  of  configurations.  There  are  no  firm  rules  at  present  and  the 
optimum  choice  is  a  strong  function  of  the  insight  of  the  particular  research 
investigator.  The  last  consideration,  proper  spin  and  symmetry  projection, 
has  proved  difficult  to  implement,  but  computer  programs  have  been  developed 
for  linear  projection  algebra  at  this  Center,  and  the  Cl  method  has  been  found 
of  demonstrable  value  in  handling  excited  states  and  dissociation  processes 
which  cannot  be  treated  with  RHF  techniques. 

Either  of  the  above  described  methods  for  ab  in it  io  calculations  reduces 
in  practice  to  a  series  of  steps,  the  most  important  of  which  are  the  evalua¬ 
tion  of  molecular  integrals,  the  construction  of  matrix  elements  of  the 
Hamiltonian,  and  the  optimization  of  molecular  orbitals  (RHF)  or  configuration 
coefficients  (Cl).  For  diatomic  molecules,  these  steps  are  all  comparable 
in  their  computing  time,  so  that  a  point  has  been  reached  where  there  is  no 
longer  any  one  bottleneck  determining  computation  speed.  In  short,  the 
integral  evaluation  involves  the  use  of  ellipsoidal  coordinates  and  the  intro¬ 
duction  of  the  Neumann  expansion  for  the  interelectronic  replusion  potential 
(Ref.  1A);  the  matrix  element  construction  depends  upon  an  analysis  of  the 
algebra  of  spin  eigenfunctions  (Ref.  15);  and  the  orbital  or  configuration 
optimization  can  be  carried  out  by  eigenvalue  techniques  (Refs.  16,  17).  All 
the  steps  have  by  now  become  relatively  standard  and  can  be  performed 
efficiently  on  a  computer  having  65,000  to  130,000  words  of  core  storage,  a 
cycle  time  in  the  microsecond  range,  and  several  hundred  thousand  words  of 
periptieral  storage. 

Both  the  RHF  and  Cl  methods  yield  electronic  wavefunctions  and  energies 
as  a  function  of  the  internuclear  separation,  the  RHF  methods  for  one  state, 
and  the  Cl  method  for  all  states  considered.  The  electronic  energies  can  be 
regarded  as  potential  curves,  from  which  may  be  deduced  equilibrium  inter¬ 
nuclear  separations,  dissociation  energies,  and  constants  describing  vibra¬ 
tional  and  rotational  motion  (including  anharmonic  and  rotation-vibration 
effects).  It  is  also  possible  to  solve  the  Schrodinger  equation  for  the  mo¬ 
tion  of  the  nuclei  subject  to  the  potential  curves,  to  obtain  vibrational 
wavefuntions  for  use  in  transition  probability  calculations.  The  electronic 
wavefunctions  themselves  can  be  used  to  estimate  dipole  moments  of  individual 
electronic  states,  transition  moments  between  different  electronic  states, 
and  other  properties.  While  all  of  the  calculations  described  in  this 


paragraph  have  been  carried  out  on  some  systems,  the  unavailability  of  good 
electronic  wavefunctions  and  potential  curves  has  limited  actual  studies  of 
most  of  these  properties  to  a  very  small  number  of  molecules. 
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METHOD  OF  APPROACH 


Central  to  these  studies  are  the  actual  quantum-mechanical  calculations 
which  must  be  carried  out  for  the  molecular  species.  For  added  clarity, 
various  aspects  of  these  calculations  are  discussed  in  individual  subsections 


1.  Levels  of  Approximation 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  near 
adequacy  of  a  minimum  Slater-type-orbital  (STO)  basis  for  constructing  molec¬ 
ular  wavef unctions  (Refs.  18  and  19).  This  means  inner-shell  and  valence- 
shell  STO's  of  quantum  numbers  appropriate  to  the  atoms  (Is,  2s,  2p  for  C,  N, 
0;  Is  for  H) .  The  main  deficiency  of  the  minimum  basis  set  is  in  its  inabil¬ 
ity  to  describe  polarization  of  :r  orbitals  in  atoms  adjacent  to  H  atoms,  and 
successful  calculations  usually  result  if  one  (or  a  set)  of  p  orbitals  is 
supplied  for  each  H  atom.  Values  of  the  screening  parameters  c,  for  each  orbi 
tal  can  either  be  set  from  atomic  studies  or  optimized  in  the  molecule;  the 
later  approach  is  indicated  for  studies  of  maximum  precision.  When  high 
chemical  accuracy  is  required,  as  for  the  detailed  studies  of  the  ground 
state  of  a  system,  a  more  extended  basis  should  be  used.  Double-zeta  plus 
polarization  functions  or  optimized  MO's  usually  are  required. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally 
flexible  manner.  This  implies  the  construction  of  Cl  wavef unctions  with  all 
kinds  of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons 
into  overall  states  can  adjust  to  an  optimum  form  at  each  geometrical  confor¬ 
mation  and  for  each  state.  Except  when  well-defined  pairings  exist  for  as 
many  electrons  as  possible,  a  single-configuration  study  (even  of  Hartree- 
Fock  quality)  will  be  inadequate. 


2.  Spin  and  Symmetry 

Proper  electronic  states  for  systems  composed  of  light  atoms  should 
possess  definite  eigenvalues  of  the  spin  operator  S^  as  well  as  an  appropriat 
geometrical  symmetry.  The  geometrical  symmetry  can  be  controlled  by  the 
assignment  of  orbitals  to  each  configuration,  but  the  spin  state  must  be  ob¬ 
tained  by  a  constructive  or  projective  technique.  Formulas  have  been 
developed  (Ref.  15)  for  projective  construction  of  spin  states  from  ortho¬ 
gonal  orbitals,  and  programs  implementing  these  formulas  have  for  several 
years  been  in  routine  use  at  UTRC. 
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One  of  the  least  widely  appreciated  aspects  of  the  spin-projection 
problem  is  that  the  same  set  of  occupied  spatial  orbitals  can  sometimes  be 
coupled  to  give  more  than  one  overall  state  of  given  S  quantum  number.  It  is 
necessary  to  include  in  calculations  all  such  spin  couplings,  as  the  optimum 
coupling  will  continuously  change  with  changes  in  the  molecular  conformation. 
This  is  especially  important  in  describing  degenerate  or  near-degenerate 
excited  electronic  states. 


3.  Method  of  Ab  Initio  Calculation 

A  spin-free,  nonrelativistic,  electrostatic  Hamiltonian  is  employed  in 
the  Born-Oppenheimer  approximation.  In  systems  containing  atoms  as  heavy  as 
N  or  0,  this  approximation  is  quite  good  for  low-lying  molecular  states.  For 
a  diatomic  molecule  containing  n  electrons,  the  approximation  leads  to  an 
electronic  Hamiltonian  depending  parametrically  on  the  internuclear  separa¬ 
tion,  R: 


J(  R) 


Z*Z 


A 


♦  z 


l>j>l 


(1) 


where  ZA  and  Zfi  are  the  charges  of  nuclei  A  and  B,  and  riA  is  the  separation 
of  electron  i  and  nucleus  A.  A  is  in  atomic  units  (energy  in  Hartrees,  length 
in  Bohrs) . 

Electronic  wavefunctions  ^(R)  are  made  to  be  optimum  approximations  to 
solutions,  for  a  given  R,  of  the  Schrodinger  equation 


M R)*(R)=  E  (R)*(R) 


by  invoking  the  variational  principle 


(2) 


„  „  /*  (RW  R)v  R)tfT 

8 W (R) =  8  - — - 

JV(R)¥(R)o'r 


(3) 


The  integrations  in  Eq.  (3)  are  over  all  electronic  coordinates  and  the  sta¬ 
tionary  values  of  W(R)  are  approximations  to  the  energies  of  states  described 
by  the  corresponding  tp ( R) .  States  of  a  particular  symmetry  are  studied  by 
restricting  the  electronic  wavefunction  to  be  a  projection  of  the  appropriate 
angular  momentum  and  spin  operators.  Excited  electronic  states  corresponding 
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to  a  particular  symmetry  are  handled  by  construction  of  configuration-interaction 
wavefunctions  of  appropriate  size  and  form. 

The  specific  form  for  sp ( R)  may  be  written 


*(R)  =  I  c„  ♦u(R)  (4^ 

H- 

where  each  ip  (R)  is  referred  to  as  a  configuration,  and  has  the  general  struc- 


where  each  iji  ^  is  a  spatial  orbital,  /}  is  the  antisymmetr iz ing  operator,  0s 
is  the  spin-projection  operator  for  spin  quantum  number  S,  and  0^  is  a  product 
of  a  and  3  one-electron  spin  functions  of  magnetic  quantum  number  M.  No  re¬ 
quirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbital,  so 
Eqs.  (4)  and  (5)  can  describe  a  completely  general  wavefunction . 

In  Hartree-Fock  calculations  ip  ( R)  is  restricted  to  a  single  t|»  which  is 
assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 
orbitals  \Jj  .  are  then  selected  to  be  the  linear  combinations  of  basis  orbi¬ 
tals  best  satisfying  Eq.  (3).  Writing 


^ui  =  £  avi  X  i 


the  a^^  are  determined  by  solving  the  matrix  Hartree-Fock  equations 


?  FX„°.MS*i£  Sx,°,,  !e0Ch  X> 


where  e<  is  the  orbital  energy  of  ip,.^ 

x  " -  MX# 


The  Fock  operator  F^v  has  been  thoroughly  discussed  in  the  literature 
(Ref.  12)  and  depends  upon  one-  and  two-electron  molecular  integrals  and  upon 
the  ay£.  This  makes  Eq.  (7)  nonlinear  and  it  is  therefore  solved  iteratively. 
UTRC  has  developed  programs  for  solving  Eq.  (7)  for  both  closed  and  open-shell 
systems,  using  basis  sets  consisting  of  Slater-type  atomic  orbitals.  Examples 
of  their  use  are  in  the  literature  (Ref.  6). 
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In  configuration  interaction  calculations,  the  summation  in  Eq.  (4)  has 
more  than  one  term,  and  the  cp  are  determined  by  imposing  Eq.  (3),  to  obtain 
the  secular  equation 


£  (Hjuv  w S/xv )  c j,  =  0  (each  fj,) 


(8) 


where 

r  / (R)  J  (R)  (R)  OfT 

Spt,=f  *£(*)*  v(R)dr  <9> 

Equation  (7)  is  solved  by  matrix  diagonalizat ion  using  either  a  modified 
Givens  method  (Ref.  16)  or  a  method  due  to  Shavitt  (Ref.  17). 

The  matrix  elements  Hpv  and  Spv  may  be  reduced  by  appropriate  operator 
algebra  to  the  forms 


H 


-l 


0  s  P 


V  (Zi.R) 


J  (R)P 


n  (l.r!) 

i=i  / 


(10) 


0m  U.p  OOV  (~',r) 


(ii) 


where  P  is  a  permutation  and  ep  its  parity.  The  sum  is  over  all  permutations 
<6m!  sPl0M>  1S  a  Sanibel  coefficient"  and  the  remaining  factors  are  spatial 
integrals  which  can  be  factored  into  one-  and  two-electron  integrals.  If  the 
are  orthonormal,  Eqs.  (10)  and  (11)  become  more  tractable  and  the  H  y  and 
S^iV  may  be  evaluated  by  explicit  methods  given  in  the  literature  (Ref.  !l3)  . 
Computer  programs  have  been  developed  for  carrying  out  this  procedure,  and 
they  have  been  used  for  problems  containing  up  to  40  total  electrons,  10  un¬ 
paired  electrons,  and  several  thousand  configurations. 


The  Cl  studies  described  above  can  be  carried  out  for  any  othronormal 
set  of  i|»  ^  for  which  the  molecular  integrals  can  be  calculated.  Programs 
developed  by  UTRC  make  specific  provision  for  the  choice  of  the  \p  as  Slater- 
type  atomic  orbitals,  as  symmetry  molecular  orbitals,  as  Hartree-^ock  orbi¬ 
tals,  or  as  more  arbitrary  combinations  of  atomic  orbitals. 
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4.  Mole.cular  Integrals 


The  one-  and  two-electron  integrals  needed  for  the  above  described 
method  of  calculation  are  evaluated  for  STO's  hv  methods  developed  by  the 
present  investigators  (Ref.  20).  All  needed  computer  programs  have  been 
developed  and  fully  tested  at  UTRC. 


5.  Configuration  Selection 

Using  a  minimum  basis  plus  polarization  set  of  one-electron  functions, 
a  typical  system  can  have  of  the  order  of  10^  configurations  in  full  Cl 
(that  resulting  from  all  possible  orbital  occupancies).  It  is  therefore 
essential  to  identify  and  use  the  configurations  describing  the  significant 
part  of  the  wavefunct ion.  There  are  several  ways  to  accomplish  this  objec¬ 
tive.  First,  one  may  screen  atomic-orbital  occupancies  to  eliminate  those 
with  excessive  formal  charge.  Alternatively,  in  a  molecular-orbital  frame¬ 
work  one  may  eliminate  configurations  with  excessive  numbers  of  anti-bonding 
orbitals.  A  third  possibility  is  to  carry  out  an  initial  screening  of  con¬ 
figurations,  rejecting  those  whose  diagonal  energies  and  interaction  matrix 
elements  do  not  satisfy  significance  criteria.  Programs  to  sort  configura¬ 
tions  on  all  the  above  criteria  are  available  at  UTRC. 

Other,  potentially  more  elegant  methods  of  configuration  choice  involve 
formal  approaches  based  on  natural-orbital  (Ref.  21)  or  multiconfiguration 
SCF  (Ref.  8)  concepts.  To  implement  the  natural-orbital  approach,  an  initial 
limited-CT  wave-function  is  transformed  to  natural-orbital  form,  and  the 
resulting  natural  orbitals  are  used  to  form  a  new  Cl.  The  hoped-for  result 
is  a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller  number 
of  significant  terms.  The  multiconfiguration  SCF  approach  is  more  cumbersome, 
but  in  principle  more  effective.  It  yields  the  optimum  orbital  choice  for  a 
pre-selected  set  of  configurations.  This  approach  works  well  when  a  small 
number  of  dominant  configurations  can  be  readily  identified. 

It  should  be  emphasized  that  the  problem  of  configuration  choice  is  not 
trivial,  and  represents  an  area  of  detailed  study  in  this  research.  The 
existence  of  this  problem  causes  integral  evaluation  to  be  far  from  a  unique 
limiting  factor  in  the  work. 


6.  Charge-Transfer  Calculations 

Even  though  low-energy  atom-atom  reactions  play  an  important  role  in  many 
physical  processes,  until  recently  comparatively  little  effort  had  been  de¬ 
voted  to  acquiring  a  knowledge  of  the  appropriate  cross  sections.  In  the 
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past,  both  theoreticians  and  experimental  physicists  have  found  it  easier  to 
study  high-energy  collisions.  At  these  energies,  the  two  colliding  particles 
preserve  their  identities,  and  it  is  possible  to  treat  the  interaction  between 
them  as  a  perturbation.  There  is  no  guarantee  that  this  procedure,  known  as 
the  Born  approximation,  will  always  converge  to  the  correct  result  (Ref.  22). 

As  the  energy  of  the  colliding  particles  decreases,  it  is  necessary  to  take 
account  of  the  distortion  these  particles  undergo  during  the  collision.  The 
method  of  perturbed  stationary  states  was  developed  for  calculating  charge 
transfer  and  electronic  excitation  cross  sections  in  relatively  slow  collisions 
between  heavy  particles  (Ref.  23).  The  method  has  been  presented  in  both 
wave  and  impact  parameter  formalisms.  In  the  first  of  these,  the  entire  sys¬ 
tem  is  treated  quantum-mechanically ,  while  in  the  latter  the  nuclei  are 
assumed  to  behave  as  classical  particles,  traveling  along  straight  line  tra¬ 
jectories,  and  the  time-dependent  Schrodinger  equation  is  solved  to  calculate 
the  probability  of  various  types  of  electronic  transitions  (Refs.  24,  25  and 
26).  Forcing  the  particles  to  travel  along  straight  lines  limits  the  validity 
of  the  impact  parameter  method  to  collisions  of  several  hundred  electron-volts 
or  greater  (Ref.  27).  The  wave  formulation  of  the  method  of  perturbed  sta¬ 
tionary  states  appears  to  be  one  practical  method  of  calculating  thermal- 
energy  charge-exchange  cross  sections.  A  semi-classical  close-coupling  met¬ 
hod  (Ref.  28),  based  on  an  average  scaling  procedure,  also  offers  utility  for 
low  to  intermediate  collision  energies. 

The  use  of  the  wave  formalism  to  study  charge-transfer  reactions  at 
thermal  energies  dates  back  to  Massey  and  Smith's  pioneering  study  of  He+  + 

He  thermal  low-energy  scattering  (Ref.  29).  Strictly  speaking,  their  theory 
is  applicable  only  to  resonant  charge  transfer  reactions,  A+A+  -*•  A  +A; 
however,  it  can  be  generalized  to  study  nonresonant  charge-transfer  reactions 
as  well  as  charge  transfer  into  excited  states.  If  the  origin  of  coordinates 
is  located  at  the  center  of  mass  of  the  nuclei  of  the  colliding  atoms,  the 
Schrodinger  equation  for  an  atomic  collision  can  be  written  in  center  of  mass 
coordinates  as 


(H-E)  X  (r.R) 


2  M 


2M 


I  I 

'  i 


V,  •  V,  -  E  )  X  (r.R)  :  0  (12) 


where  r  represents  the  position  of  the  electrons,  R  is  the  vector  joining 
A  to  B,  H0  is  the  Hamiltonian  for  the  system  when  the  nuclei  are  held  fixed, 

M  is  the  reduced  mass  of  the  two  nuclei,  M  is  the  sum  of  the  nuclear  masses, 
and  E  is  the  internal  energy  of  the  system,  including  the  electronic  energy. 
Ignoring  heavy  particle  kinetic  energy  terms  in  the  center  of  mass  system 
results  in  a  modified  form  of  the  adiabatic  approximation  (Ref.  30)  and  yields 
perturbed  molecular  eigenfunctions  iji  (r,R)  which  satisfy  the  equation 


(H°  -  in  V,.  vpi  ff.s.  <  a 

Here  en(R)  is  an  electronic  energy  level  of  the  molecule  perturbed  somewhat 
by  the  appearance  of  the  cross  terms.  The  wave  function  describing  the 
colliding  system,  X(r,R)  can  be  expanded  as 


X  (r  ,R)  =  Z  (f,R)Fn(R) 

n  '* 


(14) 


The  various  scattering  cross  sections  are  determined  by  the  asymptotic  behavior 
of  the  Fn(R).  These  functions  are  determined  by  substituting  the  expansion 
Eq.  (14)  into  Eq.  (12).  Making  use  of  the  orthogonality  of  the  molecular 
eigenfunctions,  it  is  easy  to  derive  the  following  set  of  coupled  differential 
equations  for  the  Fn(R) 


£  t, i 

m 


(15) 


*  2  *m>  *  VR  •>,<*>} 

2  2 

where  Vn(R)  and  h  /2  M  are  the  potential  and  kinetic  energies  of  particles 
in  the  n1-^1  channel.  Many  of  the  difficulties  associated  with  trying  to  cal¬ 
culate  thermal  energy  cross  sections  emanate  from  trying  to  derive  and  solve 
this  infinite  set  of  coupled  partial  differential  equations.  Until  recently, 
the  biggest  obstacle  in  the  calculation  of  low-energy  cross  sections  has  been 
the  inability  of  theorists  to  develop  accurate  molecular  eigenfunctions.  For 
those  problems  for  which  the  molecular  eigenfunctions  were  available,  the 
agreement  between  theory  and  experiment  has  been  very  reasonable  for  the 
amount  of  computational  effort  involved.  This  is  true  for  spin  exhange  and 
excitation  transfer  reactions  as  well  as  resonance  charge  transfer  reactions 
(Refs.  31-35). 


Although  the  molecular  wave  functions  and  energies  available  in  the  past 
have  not  been  sufficiently  accurate  to  permit  extension  of  the  wave  formalism 
to  systems  having  more  than  about  four  electrons,  recent  advances  in  calcula- 
tional  techniques,  especially  for  two-center  systems,  have  largely  overcome 
this  problem.  In  particular,  recent  studies  have  demonstrated  the  possibility 
of  producing  highly  accurate  adiabatic  electronic  wave  functions  for  systems 
containing  as  many  as  40  electrons  (Refs.  6,  36-38).  The  availability  of 
these  small  but  flexible  wave  functions,  which  have  the  property  of  connecting 
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with  the  correct  separated  atomic  states,  increases  substantially  the 
chances  for  successful  and  pract ical  calculat ions  of  cross  sections  using  the 
perturbed  stationary  state  technique.  This  is  especially  true  of  charge 
transfer  into  excited  states,  where  a  knowledge  of  a  number  of  the  low-lying 
excited  states  is  required. 

Even  with  a  reasonable  number  of  the  molecular  eigenfunctions,  the  prob¬ 
lem  of  calculating  cross  sections  is  far  from  solved;  this  is  especially  true 
of  charge  or  excitation  transfer  into  excited  states.  Many  of  the  existing 
studies  of  symmetric  resonance  reactions  are  based  on  the  two-state  approxi¬ 
mation  which  limits  the  number  of  terms  in  the  expansion  Eq.  (14)  to  two. 
Under  these  conditions,  the  coupling  terms  on  the  right-hand  side  of  Eq.  (15) 
vanish  (Ref.  39).  We  have,  then,  only  to  solve  two  partial  differential  equa 
tions,  instead  of  a  system  of  coupled  equations.  The  situation  is  not  as 
simple,  however,  for  nonsymmetric  reactions  or  for  high-energy  collisions 
involving  excitation.  For  these  problems,  the  coupling  terms  are  the  source 
of  the  transition  and  the  coupled  differential  equations  have  to  be  solved 
directly. 

While  a  great  deal  of  effort  has  gone  into  deriving  formal  theories  of 
inelastic  and  rearrangement  collisions,  relatively  little  work  has  gone 
into  trying  to  solve  the  resulting  equations.  This  is  especially  true  of 
low-energy  collisions,  where  the  lack  of  good  molecular  wave  functions  has 
prevented  people  from  evaluating  the  terms  coupling  the  different  channels. 
Previous  studies  which  were  part  of  the  UTRC  research  program  in  the  elec¬ 
tronic  structure  of  atoms  and  molecules  have  been  devoted  to  calculating 
matrix  elements  similar  to  some  of  the  terms  coupling  the  electronic  and 
nuclear  motion.  Some  of  the  required  work  involves  calculating  derivatives 
of  the  electronic  wave  functions  with  respect  to  the  internuclear  distance. 
This  task  is  made  simpler  by  the  use  of  compact  but  flexible  wave  functions 
such  as  those  studied  previously  at  this  center. 

The  biggest  obstacle  to  calculating  low-energy  cross  sections  is  the 
solution  of  the  infinite  system  of  coupled  partial  differential  equations 
describing  the  scattering  (Eq.  15).  The  physics  of  the  problem  usually 
serves  as  a  guide  to  truncating  these  equations  to  a  system  of  finite  order. 

A  partial  wave  expansion  of  the  F  (R)  leaves  a  large  number  of  sets  of 
coupled  differential  equations  (Refs.  40  and  41).  Since  the  number  of  equa¬ 
tions  increases  as  the  collision  energy  increases,  there  is  no  single  method 
for  solving  these  equations.  At  thermal  energies,  a  direct  numerical  integra 
tion  of  these  equations  is  feasible  (Refs.  42  and  43).  At  higher  energies, 
when  inelastic  collisions  and  charge  transfer  into  excited  states  becomes 
important,  the  trajectories  of  the  Incident  and  scattered  particles  may  bo 
nearly  classical  (Ref.  44).  Under  these  conditions,  it  is  often  possible  to 
use  the  W.  K.  B.  wave  function  to  obtain  approximate  solutions  to  these  equa¬ 
tions  (Ref.  45). 
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7.  Rotationally  Induced  Transition 


N+  +  0 


the 


which  are  important  in  the 


Our  studies  of  the  electronic  states  of  NO 
4- 

+  N  +^0  charge  transfer  reaction  led  primarily  to  consideration  of 
l  and  n  states  of  this  molecular  ion.  The  dominant  reactant  channel 
was  found  to  be  the  lowest  ^fi  state, ^which  is  bound,  and  the  connecting 
charge  transfer  channel  which  has  1  symmetry.  For  near-adiabatic  colli¬ 
sions,  this  type  of  system  can  be  handled  using  a  semi-classical  approach, 
directly  solving  the  second  order  differential  equations  which  couple  the 
adiabatic  levels,  here  of  different  angular  momentum. 


The  adiabatic  states  R)  and  ^(rjR)  are  eigenstates  of  the  non- 
relativistic  Hamiltonian,  ^  ,  for  a  fixed  internuclear  separation,  R.  We 
have 


U,R)  1  =  1,2 


(16) 


In  the  two-state  approximation  we  have,  using  impact  parameter  formulation. 


♦  =  C,(t)  (,r,R) 

+  c2(t)  *2  (r  ,r) 


-•A  Jq  ( |dr 

8 

-<A  /o’  *2dT 


(17) 


where  R  =  R(t)  and  is  defined  by  the  collision  trajectory. 

We  now  require  the  inner  product  of  i|>  and  ^  to  vanish  over  the  electron¬ 
ic  coordinates 


♦1  ^  ~  itl  I  *  >  =0 


at 


<4>z  cM  -  lt>  -L-  ^  >  rO 


at 


(18a) 

(18b) 


Combining  Eqs.  (16),  (17),  and  (18),  we  have 


?C,  . 

TT 


d<t>,  1  -M  *idT 


at 


>  8 


(19a) 


+  C*  I  “if 


a<^>2  ^-i/Vi/o  -  0 
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where  Vg  is  the  angular  velocity  of  the  internuclear  axis  and  LT  is  the  angu¬ 
lar  momentum  coupling  operator.  The  first  term  on  the  RHS  of  Eq.  (20)  leads 
to  the  well  known  Landau-Zener  (Ref.  46)  solution  for  states  of  identical 
molecular  symmetry.  For  such  cases  the  second  term  in  Eq.  (20)  vanishes. 


For  the  system  under  study  here,  only  the  second  term  on  the  RHS  of 
Eq.  (20)  contributes  in  Eqs.  (19)  and  we  are  led  to  the  following  coupled 
equations. 


c2  <<MLT  l^g>  *  * 


(21a) 


302 

at 


i  <+J|wlv't/>~“,at  =0 


(21b) 


Assuming  a  linear  dependence  of  e  on  R  near  the  crossing  point  of  the  colli¬ 
sion,  we  have 


=bMR-Rn!  =  at 


(22) 


a  :  X&  d(*2  ~C|)  (23) 

"  dR  Rx 

Assuming  also  that  <  | 1 4> 2  *  “  <LT>  is  essentially  constant  over  the 

dominant  region  of  the  collision,  we  have 
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(24a) 


i£i  „  !!»m>  c  /i/o'(‘2-,')dT;0 
at  ti 


l£l  +  iw<lt>  c  / 1. 4 

at  t>  ' 


«2)dr 

=0 


(24b) 


Equations  (24a)  and  (24b)  have  been  solved  numerically  by  Russek  (Ref. 
47)  for  the  case  of  curve-crossing  states  with  AA  =  ±  1.  For  small  values  of 
angular  velocity,  Russek  shows  that  Eq.  (24)  reduces  to  a  standard  Landau— 
Zener  form.  For  large  velocities,  the  general  solution  of  Eq.  (24)  must  be 
employed . 
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8.  Dissociative-Recombination  Rate 


The  theory  of  the  capture  of  an  electron  by  a  positive  molecular  ion  has 
been  developed  for  both  direct  attachment  processes  (Ref-  48),  and  for  several 
possible  indirect  processes  which  involve  the  formation  of  an  intermediate 
excited  Rydberg  state  belonging  to  the  molecular  ion  core  configuration 
(Refs.  49-50).  Let  (r.  Re)  represent  the  continuum  wavefunction  describ- 

'V 

ing  the  free  electron  with  energy  e  plus  the  n-1  electron  bound  system  of 
the  molecular  ion.  Let  (r,  R)  represent  the  wavefunction  for  an  eigen¬ 
state  of  the  neutral  molecule  which  can  be  written,  in  the  Born-Oppenheimer 
approximation,  as  the  product  of  an  electronic  and  nuclear  wavefunction  in 
the  form 


*d(r.  R)  ^ela.R)<d(R) 


(25) 


The  cross  section  for  electron  capture  in  the  dissociative-recombina¬ 
tion  is  determined  by  the  asymptotic  form  of  the  nuclear  wavefunction  £;d(R). 
Let  the  incident  flux  of  electrons  per  unit  area  be 


(26) 


where  <5e  is  a  measure  of  the  energy  spread  of  the  incoming  beam  and  k  is  the 
wave  number  associated  with  the  incoming  electron.  The  outgoing  flux  of 
atoms  is  given  by 


LkS« 


(27) 


where  is  the  relative  velocity  of  the  separating  neutral  atoms.  Equations 
(26)  and  (27)  yield  for  the  cross  section 


P  rr*- V  R  | 

a  (<■  ,i'.  n - T  >  (R) 

v—  v  “k1  1 


(28) 


Equation  (28)  must  finally  be  averaged  over  all  rotational  orientations  as 


a(f): 


j_.  r 

4  7 T  J 


SI 


(29) 


Assuming  that  the  total  cross  sections  can  be  written  as  the  product  of  a 
resonant  capture  cross  section  and  a  survival  factor,  in  the  form 

<r(f) -c7caD  (e)S(f)  (30 
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we  have 


*•>«•  0  fegt)  <?, 


(31) 


where  the  factor  is  the  ratio  of  the  electronic  degeneracies  of  the 

neutral  and  ionic  states,  respectively.  The  capture  width,  rc(R,e),  is 
determined  by  integrating  over  the  coordinates  of  the  electronic  wavefunctions 


.  iK«' , 


<J Ci.dXKo  tH)\J 


I'  Cl 


/el  ( J3  ,rY>2 


(32) 


Ttius,  the  capture  cross  section  can  be  written  as 


rr 


2D 


„  .  (0*2-}  1  .<"7  +  (r<  r:/<?(Mi£  (R)>2 

1  11  1  ,r.(  UA8*J  ,yfT  .‘C  ^ 


(33) 


or  equivalently 

a  (e' 


12  <34) 

The  nuclear  wavefunction,  is  energy  normalized  such  that  asymptotically 


*AB 


2  H- 


1/2 


x  ~T  )  Sin (kr  +  n) 

f\  rrk  ' 


(35) 


Equation  (34)  can  be  cast  into  the  computational  form 
crU.V)  cm3/sec 


.38188  x  10 

f(O.u) 


-is 


jAB 


uj 


’A8* 


(36) 


where  the  electron  energy,  r,  and  the  capture  width  rc  are  in  atomic  units. 

Assuming  a  Maxwellian  temperature  distribution  for  the  electrons,  the 
rate  coefficient  can  be  written  as 


■  CD 


av(Te): 


1.  1  ,  f  (T  (e,v)  X  vei  x  e 

7T  (KTe,3/l  Jo 


f/KTg  \/i 


(37) 
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or  equivalently 


a  v  (  Tp  )  - 


'rfm  (k\.) 


■v  v .  /, 


<o 


(€,v)  . 


t  >  *■ 


(U 


(38) 


The  capture  width  !’  can  be  calculated  by  examining  the  high  numbers  of  the 
Rydberg  series  of  neutral  states  which  have  the  structure  of  a  ground  state 
molecular  ion  coupled  to  an  electron  in  a  diffuse  hydrogenic  orbital  of  large 
effective  principal  quantum  numbers.  A  direct  calculation  of  the  capture 
width  involves  the  knowledge  of  the  continuum  wavef unction  for  the  electron 
as  a  function  of  the  interpar tide  coordinates.  This  approach  is  computa¬ 
tionally  very  difficult  and  at  the  present  stage  of  development  would 
probably  lead  to  errors  of  at  least  the  same  magnitude  as  an  extrapolation 
procedure  of  the  corresponding  neutral  Rydberg  states.  An  alternate  pro¬ 
cedure  is  given  here  which  involves  analyzing  only  bound  state  functions. 

Evaluation  o_f  Elec. tronicjrransi_t_ion  Mat r jLx_Elemen  t s_ 

One  of  the  most  difficult  tasks  in  the  calculation  of  the  cross  section 
for  dissociative-recombination  is  the  evaluation  of  electronic  matrix  elements 
of  the  form  R)  |  d  .  -e  U'fn+(r»  R)f(t,  r,  R)>  where  <*4  .  is  the  total 

electronic  Hamiltonian,  ^‘;^(r,  R)  is  the  electronic  wavefunction  for  the 

el  ^  % 

valence  state,  and  ti>.r,+  (r,  R)f(c,  r,  R)  is  the  continuum  wavefunction  with 
total  electronic  energy,  c.  Unlike  the  which  correspond  to  bound  elec¬ 
tronic  valence  states,  the  (r,  R)f(e,  r,  R)  do  not  vanish  for  large 

electron-molecular  ion  separations  but  have  an  oscillatory  behavior.  The 
direct  evaluation  of  these  wavef unctions  involves  the  solution  of  a  very 
complicated  scattering  problem,  and  cannot  be  directly  accomplished  by  the 
standard  variational  methods.  We  have  developed  two  independent  methods 
to  evaluate  these  matrix  elements  that  do  not  require  the  direct  solution 
of  the  electron-molecular  ion  scattering  problem. 

Since  it  is  a  straightforward  process  to  evaluate  matrix  elements  of  the 
form  j  ei-e  1 ^  where  E’e*(~En*)  is  the  wavefunction  for  a  bound 

Rydberg  state,  one  way  of  estimating  the  threshold  value  of  the  continuum 
analogue  of  these  integrals  is  to  analytically  continue  the  bound  state 
matrix  elements  into  the  continuum.  In  performing  this  analytic  continuation, 
one  has  to  be  careful  to  include  a  factor  corresponding  to  the  density  of 
states,  oCE^a)  -  dn/dE.  While  this  factor  is  implicitly  included  in  the 
normalization  constant  for  the  continuum  wavefunctions ,  it  has  to  be 
explicitly  calculated  for  ttie  bound  state  functions.  To  a  first  approxi¬ 
mation,  the  energy  of  the  Rydberg  states  is  proportional  to  therefore, 

the  density  of  states  increases,  approximately,  as  a ^ .  On  tiie  basis  of 
these  approximations,  we  have  been  able  to  estimate  the  bound-continuum 
matrix  elements  for  low  scattering  energies. 
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We  have  derived  a  second  extrapolation  procedure  to  calculate  the 
bound- continuum  electronic  matrix  elements  at  higher  energies.  A  common 
procedure  in  theoretical  physics  for  converting  a  continuum  or  scattering 
problem  to  a  finite  eigenvalue  problem  is  to  enclose  the  system  in  a 
spherical  shell  of  finite  volume  and  require  the  wavefunction  to  vanish  at 
the  boundaries.  Depending  on  the  form  of  the  scattering  potential,  there 
will  be  a  finite  number  of  bound  or  negative  energy  states,  and  an  infinite 
number  of  positive  energy  states,  all  satisfying  the  same  boundary  conditions. 
If  the  radius  of  the  spherical  shell  is  large  enough,  these  bound  state  wave- 
functions  will  be  a  reasonable  approximation  to  the  exact  wavefunctions  for 
the  bound  states.  For  any  bound  state  and  the  infinite  number  of  discrete 
states,  it  is  possible  to  estimate  the  value  of  any  bound-continuum  matrix 
elements  by  plotting  values  of  the  product  of  the  square  of  these  matrix 
elements,  times  a  discrete  density  of  states  function  p(E([*)  =  j  1/ (E  j  . 

If  the  radius  of  our  imaginary  sphere  is  large  enough,  the  curve  generated 
from  this  procedure  should  give  reasonable  values  for  these  matrix  elements. 

We  have  derived  a  many  electron  version  of  this  theory.  A  finite  set 
of  bound  state  functions,  which  has  a  large  continuum  contribution,  is 
diagonalized  with  respect  to  the  operator ,  Jfe\ .  This  special  set  of 
Rydberg  functions  is  then  used  to  generate  a  series  of  curves  for  values  of 
the  matrix  elements  <\bv  Wei_E:  A  set  of  estimated  values  of  these 

matrix  elements  is  given  in  Table  A  for  s-wave  (sigma)  scattering  states  of 
NO. 


Energy* 

.73 

2.39 

4.16 

6.32 

10.91 

*energy  in  atomic  units  relative 
Indirect  Models 


TABLE  A 

?(E)*  <^t2IqIT)Me|-ekRyd>2 

8.31  E-06 
1.43  E-06 
2.36  E-06 
2.72  E-06 
6.52  E-08 


An  indirect  model  for  calculating  electron-molecular  ion  recombination 
coefficients  can  also  be  formulated.  Here  we  postulate  the  formation  of  a 
collision  complex  (AB*)  during  the  collision  and  its  subsequent  decay  to  form 
A*,  B**,  or  AB+  +  e  (Refs.  49-51).  The  approaches  differ  by  the  nature  of 
the  (AB*)  complex.  The  high-energy  approximation  assumes  that  all  of  the 
energy  of  the  incident  electron  goes  into  exciting  one  of  the  electrons  in 
the  AB+  core,  resulting  in  the  capture  of  the  incident  electron  into  a 
doubly  excited  state  of  AB.  This  mechanism  is  unlikely  to  be  applicable  to 


20 


AB+  recombination  at  thermal  energies  because  of  the  energy  required  to 
promote  an  electron  out  of  the  AB+  core.  A  second  approach  is  more 
applicable  to  thermal  energy  collisions,  since  it  requires  that  the  electron- 
molecular  ion  collision  results  in  a  vibrat ionally  excited  molecular  ion 
which  captures  an  electron  in  the  Coulomb  field  of  AB+.  This  collision 
complex  is  equivalent  to  a  vibrationally  or  rotationally  excited  Rydberg 
state  of  AB.  Experimental  evidence  for  the  existence  of  such  states  and 
their  autoionization  has  been  independently  obtained  from  photoionization 
spectra  of  AB,  within  a  few  electron  volts  of  threshold  (Ref.  52). 

With  this  model,  the  prediction  of  the  atomic  products  of  the  recombination 
reaction  reduces  to  calculating  the  Rydberg  states  of  AB  that  are  most  likely 
to  be  populated  during  the  initial  recombination,  and  also  those  atomic 
states  to  which  they  subsequently  decay. 

As  contrasted  with  the  direct  recombination  process  which  is  governed 
by  tile  configuration  interaction  strength,  the  nonadiabatic  coupling  of  the 
electron  and  nuclear  motion  is  the  dominant  mechanism  in  both  stages  of  the 
indirect  recombination  process.  In  the  first  stage,  it  is  the  vibrational 
(or  rotational)  excitation  of  the  AB+  core  whicli  results  in  the  initial 
capture,  of  the  incident  electron.  The  subsequent  decay  of  the  collision 
complex  is  due  to  a  transition  to  a  repulsive  molecular  state  which  crosses 
the  Rydberg  state.  It  is  then  the  nonadiabatic  coupling  of  the  electronic 
and  nuclear  motion  which  results  in  the  transition  between  these  two  molecular 
states.  Actually,  it  is  the  same  term  in  the  Hamiltonian  for  the  entire 
molecular  system  that  causes  the  key  transition  in  both  stages  of  the  recom¬ 
bination  process.  This  becomes  obvious  in  the  quantum  mechanical  formulation 
of  the  problem.  A  possible  competing  process  for  the  decay  of  the  AB*  complex 
is  the  radiative  transition  to  lower-lying  states  of  the  molecule.  The 
importance  of  this  effect  can  be  estimated  from  calculations  of  the  band 
transition  probabilities. 

Implicit  in  this  model  of  dissociative-recombination  processes  is  the 
existence  of  three  excited  electronic  states  of  AB.  Since  this  recombination 
process  is  envisioned  as  occurring  in  two  stages,  we  are  concerned  with  only 
two  of  these  states  at  a  time.  The  state  common  to  both  stages  of  the 
recombination  process  is  the  highly  excited  Rydberg  state  (AB*) ,  which  is 
bound  with  respect  to  dissociation  to  AB+  and  a  free  electron.  For  the  calcu¬ 
lation  of  the  initial  rate  of  recombination,  a  wavefunction  for  a  continuum 
Rydberg  state  is  needed.  This  wavefunction  has  the  asymptotic  behavior  of 
a  free  electron  moving  in  the  field  of  AB  . 


To  determine  the  transition  rates  in  the  two  stages  of  the  recombination 
process,  it  is  necessary  to  identify  the  terms  in  the  Hamiltonian  that  are 
responsible  for  the  transition.  The  Hamiltonian  for  the  complete  molecular 
system  with  respect  to  the  center  of  mass  is 


21 


where  R  Is  the  vector  joining  nucleus  B ,  is  the  Hamiltonian  for  the 
system  where  the  nuclei  are  held  fixed,  u  is  the  reduced  mass  of  the  two 
nuclei,  M  is  the  sum  of  the  nuclear  masses,  and  the  summations  extend  over 
all  of  the  electrons.  Since  all  of  the  electronic  eigenfunctions  under 
discussion  are  implicitly^dependent  on  the  internuclear  separation  jj^|,  it 
is  the  term  Jf  =  -(h2/2p)V^  which  plays  a  dominant  role  in  coupling  the 

electronic  and  nuclear  motion  giving  rise  to  the  appropriate  transition 
(Refs.  53-56).  This  becomes  more  obvious  when  we  write  the  total  wavefunction 
for  the  collision  complex,  AB*,  as 

^c'r  XcUe-M)  FcW 

c  c  ~  j~r  (40) 

where  ^c(*e»lRl)  is  the  electronic  wavefunction  for  the  appropriate  Rydberg 
state  of  at  internuclear  separation  | R |  and  Fc(£)  is  the  wavefunction 
associated  with  Che  vibration  and  rotation  of  the  nuclei.  In  a  similar 
manner,  the  total  wavefunction  associated  with  the  repulsive  state  is 


\  =X*(£e,|R|)F*  (R) 


where  F.(R)  is  the  wavefunction  f-  r  the  nitrogen  and  oxygen  atoms  as  they 
separate.  With  the  use  of  Fermi's  golden  rule,  the  transition  rate  for  the 
decay  of  the  collision  can  be  obtained  by  evaluating  the  matrix  elements. 


■|<kK!*>; 


The  initial  rate  of  recombination  into  the  highly  excited  Rydberg  state  of 
nitrogen  can  be  calculated  in  a  similar  manner. 

We  now  consider  the  problem  of  the  evaluation  of  the  matrix  elements 
appearing  in  Eq.  (42).  The  most  obvious  way  of  handling  this  is  to  employ 
a  direct  approach  and  evaluate  all  matrix  elements  directly  by  quadrature. 
This  reduces  the  problem  to  the  calculation  of  matrix  elements  of  the  form 
(Ref.  57). 


<  Fc  <5'  | Vr  1 X  ,  ( r,e  ,|a|) (8> 


While  this  task  is  possible,  some  of  the  physics  of  the  problem  become  obscure 
with  this  direct  approach.  This  is  more  obvious  if  Eq .  (43)  is  simplified  to 
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<x  c  (JL e  - 1 R  |)  Fc  (R )  I  V vt  x  *  ( J>  R  | )  +  2Vr  x  *  VR  F*  +  *  t  v|  F*  >  ( 44 ) 


Since  the  derivatives  of  F.(R)  with  respect  to  the  nuclear  coordinates  are 

*■  <~\j 

much  greater  than  those  of  Xc(re,R),  Eq .  (44)  reduces  to 


2  Oc  We ,  I  a|)  *4>  <Fc  (R)  VR  (R)> 


which  can  be  expanded  as 


ir>  +  <^(R)Fcfe> 


jFJ_  X  (46) 


where 


. '  2<Cxc <ie,|B|)  Tfi*  <ie' R|l> 


^,’2<Xe(ie.l5|)XlIe.|«[)> 


(47a) 


(47b) 


If  XQ  and  Xr  are  both  eigenfunctions  of  the  electronic  Hamiltonian,  i.e.. 


Xciit’  3  )^C(R)XC 


Re/  X  n,  3 


.(|r|)Xi 


(48a) 


(48b) 


Nct  (R) -  2 - -  <  v  (*3e[Rl )  — ® —  x n  X  Z 

«*(R)-«C(R)  *c  1  1  Rr}  * 


(49a) 


(49b) 
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2 

where  the  many-electron  operator  Z^/Rg  refers  to  a  coordinate  system 
centered  on  the  nucleus  B.  for  o- a  transitions,  Nc  ^  =  0  and  for  a -if 
transitions  M  =0.  The  computational  effort  is  thus  directed  to  the 
evaluation  of  the  matrix  elements  appearing  in  Eqs.  (49).  Upon  expansion 
of  the  many-electron  operators,  such  as  Zg/Rg  >  with  respect  to  the  center 
of  mass  of  the  nuclei,  the  long-range  interactions  that  appear  in  the  more 
qualitative  theories  of  dissociative-recombination  and  autoionization 
(Refs.  54  and  55)  are  recovered. 
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DISCUSSION  Of  RESULTS 


The  theoretical  research  conducted  under  this  program  was  concerned  with 
the  calculation  of  the  dissociat ive-recombinat ion  reaction  of  electron 
attachment  to  the  nitric  oxide  molecular  ion  and,  in  particular,  to  the 
determination  of  the  dissociation  products  of  this  reaction.  An  important  part 
of  this  study  was  the  determination  of  accurate  wavef unctions  and  potential 
energy  curves  for  selected  electronic  states  of  the  nitric  oxide  molecule 
and  the  evaluation  of  product  state  branching  using  a  quantum  mechanical  close¬ 
coupling  collision  code.  A  second  area  of  research  involved  a  preliminary 
investigation  of  the  most  probable  reaction  pathway  for  the  ion-molecule 
reaction  0+  +  N2  N0+  +  N. 

Dissociative-Recombination  in  N0+ 


A  theoretical  research  program  was  undertaken  and  directed  toward  the 
study  of  the  kinetics  of  low-energy,  electron-molecular,  ion  collision 
processes.  This  program  was  primarily  concerned  with  the  theoretical 
determination  of  the  branching  ratio  for  e  +  N0+  ->  N[^D,  ^S]  +  0  [^P]. 

The  first  step  in  this  research  program  was  the  accurate  determination  of 
the  electronic  states  of  NO  which  could  represent  possible  dissocation 
channels  for  the  reaction.  Symmetry  considerations  dictate  that  the  electronic 
states,  ^I+,  2H ,  are  the  most  important  channels.  States  of  negative 

parity  or  of  quartet  multiplicity  cannot  connect  with  the  initial  state  of 
e  +  N0+  in  a  first  order  approximation.  Only  spin-orbit  coupling  or  two- 
electron  excitation  matrix  elements  (both  of  which  are  weak  for  this  system) 
can  give  rise  to  second-order  coupling.  Higher  electronic  states  have  non¬ 
connecting  matrix  elements  for  the  interaction  potential.  The  method  of 
calculation  which  was  employed  for  these  electronic  states  of  fiO  (full  Cl) 
was  capable  of  yielding  wavef unctions  which  asymptotically  connect  to  the 
correct  ato  ic  limits,  but  which  indicated  a  smoothly  varying  end  increasing 
error  for  shorter  internuclear  separations  owing  to  nonoptimization  of  the 
inner-shell  basis.  Tills  error  could  be  analyzed  quantitatively  for  any  given 
molecular  symmetry  by  comparing  the  calculated  potential  curves  with  the 
known  experimental  levels  for  the  lowest  state  of  a  given  symmetry.  Application 
of  these  correction  curves  to  the  calculated  potential  curves  yields  a 
rationalized  estimate  of  the  true  behavior  of  the  excited  electronic  states 
of  any  given  symmetry.  This  procedure  is  based  on  the  assumption  that  the  full 
d  calculations  yield  a  correct  relative  picture  of  the  excited  electronic 
states.  This  implies  that  most  of  the  valence  shell  correlation  energy  has 
been  adequately  described  by  these  full  Cl  calculations. 
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The  results  of  this  study  of  the  electronic  states  of  NO  indicate  that 
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states  of  }  ,  *11 ,  and  *~A  symmetry  will  be  effective  in  low-energy  dissociative- 

recombination  of  a  +  NO  .  Calculations  were  performed  for  electronic  states  of 
NO  in  various  approximations  including  minimum  basis,  double  zeta  and  d-type 
polarization  functions.  The  final  adjusted  energies  for  2Z+,  2il ,  and  2A 
symmetries  of  NO  are  shown  in  Fig.  1  which  represents  the  most  probable  set 
of  potential  curves  which  lead  to  dissociative-recombination. 

The  dissociative-recombinations  cross-sections  are  evaluated  emploving  the 
direct  attachment  process.  This  direct  process  is  governed  by  the  matrix 
element  <^0+  *  f  |  V(R)  |  where  refers  to  the  dissociating  neutral 

molecular  state  of  NO.  Previous  calculations  have  shown  that  the  direct 
attachment  process  will  dominate  if  the  Franck-Condon  overlap  between  N0+ 
and  the  dissociating  channels  of  NO  is  favorable.  Figure  1  indicates  this 
to  clearly  be  the  case  for  this  system. 

The  calculated  rate  constants  as  a  function  of  the  electron  temperature 
for  the  dissociating  channels  2Z+,  B2n  n,  2T1  III,  2  i, ,  and  ~t>  are  given  in 
Tables  1  and  2.  These  channels  are  the  only  significant  dissociation  channels 
up  to  about  20,000°K.  These  calculated  rates  show  a  complicated  behavior  with 
vibrational  quantum  number  of  the  molecular  ion  and  tend  to  reinforce  our 
previous  conclusions  (Ref.  58)  that  simple  models  which  assume  that  recombination 
into  higher  vibrational  states  of  the  ion  can  be  neglected  are  probably 
invalid.  Here  we  find  recombination  into  the  second  vibrational  state  occurs 
at  a  larger  rate  than  into  the  ground  state.  The  more  basic  parameters  for 
this  reaction  are  the  collision  cross-sections  as  a  function  of  electron 
energy  of  the  N0+  ion.  These  are  presented  in  Table  3  for  recombination  into 
the  ground  vibrational  state. 

The  summed  contribution  of  all  dissociating  channels  is  shown  in  Tables  1 
and  2.  Only  the  and  If  states  appear  to  contribute  significantly  to  the 

overall  recombination  rate.  The  contribution  of  '  is  f>0  percent  with  the 
remainder  being  primarily  B2n  IT,  at  least  for  temperatures  below  n.  10 ,000°K. 

2 

At  higher  energies  the  n  III  becomes  increasingly  more  important  and,  in 
fact,  is  the  dominant  dissociating  channel  above  'v  2.0  eV.  This  situation 
leads  to  a  total  recombination  cross-section  with  some  indicated  structure. 

A  comparison  of  these  calculated  data  with  the  most  recent  experimental  data 
is  shown  in  Fig.  2. 

All  of  the  primary  dissociating  channels  lead  to  N  2l)  +  0  atomic  states 

')  + 

as  can  be  seen  in  Fig.  1.  The  only  ambiguity  is  for  ^T.  symmetry  where  an 
avoided  crossing  occurs  between  states  I  2I+  ]  and  A  2r.+  II.  The  former  leads 
to  N  ^S+  0  and  the  latter  to  N  2D+  0  ^P  if  diabatic  behavior  is  followed 
through  the  crossing  point.  A  bandau-Zener  analysis  at  this  crossing  point 
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indicates  that  the  diabatic  curves  should  be  followed,  but  the  calculation  is 
very  sensitive  to  the  choice  of  the  interaction  energy,  H^..  Thus  a  full 
quantum  mechanical  close-coupling  calculation  of  the  branching  was  required 
for  this  system. 

There  are  several  approaches  to  the  quantum  theory  of  reactive  scattering. 
Most  of  these  approaches  differ  in  their  treatment  of  the  coordinate  framework 
for  describing  the  collision  process  and  in  the  detailed  asymptotic  description 
of  reactants  and  products. 


In  the  present  study,  we  utilize  the  general  approach  of  Diestler  (Ref.  59) 
as  implemented  by  Schmalz,  Stechel,  and  Light  (Ref.  60)  for  the  general 
problem  of  charge  transfer  and  exchange  reactions.  In  this  work  a  non- 
orthogonal  or  diabatic  basis  set  is  utilized  to  represent  the  reactant  and 
product  states.  The  solution  of  the  coupled  Schrbdinger  equations  describing 
the  nuclear  motion  is  carried  out  numerically  using  modified  R-matrix 
propagator  techniques. 

The  reaction  under  study  can  be  written  as: 


N  [2D]  +  0  [3P]  +  NO  [A  2E+]  (50) 

NO  [A  2I+]  -*■  NO  [  I  2E+] 

NO  [I  2E+  ]  •+  N  [4S]  +  0  [3P] 

k2  ^  0.3844  eV 


The  branching  ratio  analysis  is  illustrated  in  Pig.  3  which  indicates  that  for 
dissociative-recombination  of  e  +  N0+  into  the  A  E+  channel,  we  need  to 
consider  the  equivalent  scattering  problem  shown  in  Eq .  (50)  for  scattering 
energies  5^  0.3844  eV,  a  value  corresponding  to  zero  kinetic  energy  attachment 
of  an  electron  to  N0+. 


2  + 

The  adiabatic  potential  curves  for  the  lowest  two  I  states  of  NO  were 

converted  to  diabatic  form  by  interpolation  through  the  crossing  region.  Our 

2  4- 

calculated  diabatic  potential  curves  for  E  symmetry  are  given  in  Table  4 
where  we  see  that  the  maximum  Hamiltonian  interaction  term,  j ,  is  270  cm-3 
at  an  internuclear  separation  of  2.842  bohrs.  The  detailed  shape  of  this 
diabatic  interaction  potential  is  shown  in  Pig.  4.  The  smallness  of  this 
interaction  term  indicates  that  coupling  is  weak  between  these  lowest  two  E+ 
states  and  that  diabatic  collision  behavior  is  indicated. 
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Detailed  calculations  of  the  cross-sections  for  the  reaction  illustrated 
in  Eq.  (50)  were  carried  out  for  collision  energies  0.3844  <_  E(eV)  <_  2.00. 


This  corresponds  to  electron  attachment  in  the  energy  region  0.  E(eV)  <_  1.6  eV. 

Our  calculated  collision  cross-sections  are  shown  in  Table  5  where  we  see  that 

2  2  + 

the  coupling  cross-section  between  I  zI+and  A  E+of  NO  never  exceeds 
-5  2  -22  2 

o  'v  1  x  10~J  a  (3  x  10  cm'  )  .  We  are  thus  led  to  the  prediction  that  the 

O  ?  4- 

recombination  along  the  A  I  channel  follows  a  diabatic  reaction  path  and  we 
predict  exclusively  that  N(^D)  +  0(^P)  are  the  atomic  products  of  dissociative- 
recombination  of  e  +  N0+. 


Comparison  With  Experimental  Data 


-0  39 

The  calculated  data  show  a  low-temperature  variation  as  T£  ‘  fpr  v  =  0 
and  Te  ‘  for  v  =  2,  a  result  which  is  not  in  accord  with  the  results  obtained 
in  a  microwave-afterglow-mass-spectrometer  by  Gunton  and  Shaw  (Ref.  61)  but 
is  in  excellent  agreement  with  the  experimental  data  of  Huang,  Biondi  and 
Johnson  (Ref.  62).  This  agreement  may  be  somewhat  fortuitous  since  Huang, 
et  al.  almost  certainly  measured  recombination  primarily  into  the  second 
vibrational  level  of  NO+.  Our  predicted  temperature  dependence  of  the  rate 
coefficient  for  v  =  2  is  in  somewhat  poorer  agreement  with  their  experiments. 


Another  feature  of  the  comparison  with  the  data  of  Huang,  et  al .  is  the 
indicated  strong  increase  in  the  recombination  rate  at  temperature  <  300°F. 
This  may  be  due  to  secondary  attachment  processes  (recombination  into  vibra- 
tionally  excited  Rydberg  states  of  NO)  which  will  become  more  important  at 
lower  electron  temperatures  (Ref.  63).  The  overall  rates  are  in  numerical 
agreement  within  our  calculation  uncertainty. 


The  recent  ion  recombination  studies  by  Walls  and  Dunn(Ref.  64)  show  a 
stronger  energy  dependence  in  the  measured  cross-sections  than  that  found 
theoretically  from  these  studies  for  electron  temperatures  below  about  0.1  eV. 
Examination  of  Fig.  5  shows  that  the  cross-section  should  approach  1/e 
behavior  at  about  0.04  eV  whereas  the  data  of  Walls  and  Dunn  clearly  exhibit 
a  different  energy  dependence  in  this  region.  The  major  problem  seems  to 
be  the  experimental  data  point  at  ^  0.04  eV.  A  recent  reexamination  of  these 
data  by  Walls  and  Dunn  (Ref.  65)  indicate  a  large  experimental  uncertainty 
at  these  very  low  electron  energies;  their  lower  uncertainty  nearly 
bracketing  our  calculated  data.  The  inferred  recombination  rates  from  the 
data  of  Walls  and  Dunn  show  a  very  strong  temperature  dependence  (T  ) 

arising  from  the  low  energy  dependence  of  their  measured  cross-section  data. 

The  agreement  between  calculated  and  measured  cross-sections  above  0.1  eV 
is  excellent  even  to  systematics  in  the  structure  whicli  arise  from  a  changing 
importance  of  the  dissociating  channels  as  a  function  of  increasing  electron 
energy.  The  flattening  of  the  cross-section  above  1.5  eV  is  the  result  of 
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competing  effects  between  the  B^TI  II  and  m  dissociating  channels.  The 
latter  is  becoming  more  important  in  this  energy  region  and  shows  a  weaker 
energy  dependence,  hence  a  leveling  off  effect  occurs  in  the  total  cross- 
section.  This  qualitative  agreement  tends  to  verify  that  we  are  considering 
the  correct  molecular  states  of  NO  as  the  dissociating  channels. 

The  data  of  Gunton  and  Shaw  (Ref.  61)  are  suspect  owing  to  dimer  ion 
effects.  The  importance  of  ion  clustering  at  low  temperatures  is  uncertain 
for  the  other  experiments. 

A  further  question  is  the  importance  of  autoionization,  which  may  be 
enhanced  at  the  lower  energies  for  this  reaction.  A  crude  calculation  of  the 
survival  factor  as  defined  by  Bardsley  (Ref.  49)  indicates  that  it  should  be 
close  to  unity.  Further  detailed  studies  are  required  for  a  more  quantitative 
analysis  of  the  role  of  autoionization  at  the  lower  electron  energies. 

Our  calculated  rate  coefficients  may  be  considered  to  be  valid  only  up 
to  about  10,000°K.  Above  this  temperature,  our  analysis  is  incomplete  owing 
to  increased  population  of  the  higher  vibrational  structure  and  the  neglect 
of  other  electronic  states.  The  comparison  with  the  high-temperature  data 
of  Stein  (Ref.  66)  and  Lin  and  Tease  (Ref.  65)  appear  to  indicate  a  much 
stronger  electron-temperature  dependence  than  that  found  in  these  calculations 
or  that  can  even  be  inferred  from  the  data  of  Walls  and  Dunn. 

0+_+_N2_Jon_Molecule  Reaction_ 

As  a  continuation  of  our  theoretical  research  on  atmospheric  chemistry, 
we  undertook  a  preliminary  examination  of  the  reaction  kinetics  and  branching 
ratios  for  the  ton-molecule  reaction  0+  +  N7  N0+  +  N.  The  particular 
reactions  which  are  most  important  in  N0+  formation  are: 


o+  [4s]  +  n2^  [x  1::g+]  -  no+  lx  V]  +  n  l4s]  (51) 

AE  =  +  1.1  eV 

0+  [4s]  +  N2*  [A  3ru+]  -*  N0+  [X  V]  +  N  [4s]  (52) 

AE  =  +  7.3  eV 


The  molecular  correlations  which  connect  reactant  to  product  states  are 
uncertain.  A  complete  N20  correlation  diagram  is  shown  in  Table  6  which 
lists  the  reactant  and  product  states  in  increasing  energy  without  symmetry 


considerations.  The  lowest  quartet  state  in  a  linear  (CmV)  0-N2+ 
figuration  corresponds  to  symmetry.  In  a  bent  (C  )  geometry. 


con- 

this 


lowest 
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configuration  belongs  to  ^A".  Hopper  (Ref.  68)  has  carried  out  SCF 
calculations  for  several  symmetries  of  N2O  and  has  suggested  a  complicated 
reaction  pathway  in  which  0  +  first  combine  in  a  linear  configuration 

and  then  form  a  bent  conformation  with  N  [^S]  finally  leaving  at  90°  relative 
to  the  original  collision.  Hopper  calculates  the  lowest  saddle  point  in  the 
collinear  symmetry  to  be  a  least  7  eV  above  the  0+  [^S]  +  N2  [^Sg+] 
asymptote.  This  is  energetically  inaccessible  for  low-energy  0+  +  NQ 
collisions  and  thus  he  favors  the  distorted  reaction  path  in  symmetry. 

Previous  studies  in  this  Center  (Ref.  69)  on  N02+  have  shown  that  SCF 
or  MCSCF  calculations  may  exhibit  artifically  high  reaction  barriers  in  CMV 
symmetry.  The  origin  of  the  problem  is  that  SCF  or  MCSCF  methods  generallv 
employ  maximum  pairing  in  their  MO  structure.  The  SCF  structure  of  ^O"1"  [^7~] 
can  be  written  as: 


la2  2ct2  Jo2  4o2  5o2  bo2  In'*  7o2  2rr2 


8o 


where  at  short  0+  -  N2  separations  we  can  identify  the  principal  AO  composition 
of  these  MO's  as  follows: 


la  ,  2 a  ,  3a  = 
40  = 
5o  = 
6a  = 
7a  = 
ga  = 
90  = 

It  = 
2tt  = 
3t  = 


Is 

2s0 

2sn 

(outer) 

2sn 

( inner) 

2  pa 

(NO) 

2poN 

(outer) 

2po* 

(NO) 

2pn 

(NO) 

2piry 

(inner) 

2pTT 

(NO) 

At  large  0  -  ^  separations,  the  MO  composition  for  the  2p  electrons  must  be: 


7q  =  2po  (N2) 
8a  =  2po  (0+) 

It  =  2pT  (N2) 

2-  =  2pn  (0+) 


Thus  both  the  a  and  tt  orbitals  are  completely  recoupled  spin-wise  in  going  from 
0+  +  N2  to  N  t  N0+.  in  Fig.  6  we  present  a  general  vector  coupling  diagram 
which  indicates  that  there  are  48  separate  spin-states  for  nine  electron 
quartet  coupling.  Reactions  (51)  and  (52)  involve  a  complete  spin-recoupling 
as  the  pairing  changes  in  going  from  reactants  to  products  which  suggests  that 
a  valence  bond  treatment  is  more  appropriate  than  SCF  or  limited  MCSCF 
expansions.  Previous  studies  on  the  much  simpler  surface  have  indicated 
the  importance  of  including  higher  spin-coupled  states  (Ref.  70). 
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Our  preliminary  studies  indicate  that  the  spia-recoupling  that  occurs 
in  going  from  reactants  to  products  is  not  adequately  represented  within 
SCF  or  MCSCF  frameworks.  A  calculation  for  of  ^0"*"  at  intermediate 
internuclear  spearations  yields  a  spin-correlation  error  of  6  eV.  This 
value  closely  corresponds  to  the  7  eV  barrier  found  in  the  MCSCF  calculations 
by  Hopper.  We  propose  to  examine  the  lowest  energy  potential  energy 
surface  of  N20+  within  an  ab-initio  Cl  framework  including  all  possible 
spin  and  spatial  projections  of  the  2p  electrons.  Based  on  our  previous 
studies  of  H-j  and  NC>2,  we  believe  that  a  linear  pathway  of  much  lower 
energy  than  that  predicted  from  SCF  calculations  will  be  found.  A  schematic 
of  our  predicted  reaction  pathway  is  illustrated  in  Fig.  7.  A  similar 
analysis  wil  be  carried  out  for  0+(^S)  +  (A  ^£u+)  >  where  both  doublet 

and  quartet  channels  will  be  examined. 
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Figure  1  Potential  Energy  Curves  of  NO  Leading  to  Dissociative  Recombination  in  NO  + 
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FIGURE  2.  DISSOCIATIVE  RECOMBINATION  RATE  COEFFICIENT,  cm3/SEC 
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Figure  3  Branching  Ratio  Analysis 
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Figure  4  Interaction  Potential  lor  NO  (A  2  £+  -|  2  £+) 


42 


RECOMBINATION  CROSS  SECTION  (cm2) 


•♦NO* O'*  01 


N  ♦  0 


ELECTRON  ENERGY  (eV) 

FIGURE  5.  TOTAL  RECOMBINATION  CROSS  SECTION,  cm2 
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Figure  7  Co-Linear  Reaction  Pathway  for  0+  +N2 
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TABLE  1 


CALCULATED  DISSOCIATIVE-RECOMBINATION  RATE  COEFFICIENTS 
a  (Te,  v  =  0) ,  cm^/sec 
e  +  N0+  (v  =  0)  ->  N  +  0 


T  °K 

a 

a  B  2n 

2 

a  n  III 

a  B'  2 A 

a  0 

a  Total 

200 

9.216-08 

4.437-08 

1.615-09 

4.423-09 

1.342-10 

1.427-07 

298.16 

7.806-08 

3.669-08 

1.366-09 

3.614-09 

1.122-10 

1.198.07 

400 

6.960-08 

3.194-08 

1.221-09 

3.114-09 

9.895-11 

1.060-07 

600 

5.969-08 

2.640-08 

1.071-09 

2.538-09 

8.420-11 

8.978-08 

800 

5.348-08 

2.310-08 

9.975-10 

2.195-09 

7.599-11 

7.985-08 

1000 

4.918-08 

2.085-08 

9.610-10 

1.962-09 

7.090-11 

7.302-08 

2000 

3.811-08 

1.526-08 

9.888-10 

1.377-09 

6.324-11 

5.580-08 

4000 

2.903-08 

1.122-08 

1.416-09 

9.256-10 

7.218-11 

4.266-08 

6000 

2.347-08 

9.167-09 

2.027-09 

6.993-10 

8.009-11 

3.544-08 

8000 

2.020-08 

7.919-09 

2.581-09 

5.502-10 

8. 11 3-11 

3.133-08 

10000 

1.715-08 

6.892-09 

2.970-09 

4.517-10 

7.792-11 

2.754-08 

TABLE  2 

CALCULATED  DISSOCIATIVE-RECOMBINATION  RATE  COEFFICIENTS 


a 

(Te,  v  =  2), 

cm3/sec 

e 

+  NO+  (v  =  2) 

-»  N  +  0 

T  °K 

a  V 

a  B2n 

a  2n  III 

a  B'  ^ A 

2 

a _ $ 

a  Total 

200 

6.461-08 

1.527-08 

1.288-07 

1.498-09 

1.137-10 

2.103-07 

298-16 

5.406-08 

1.277-08 

1.057-07 

1.194-09 

I 

8.781-11  , 

1.738-07 

400 

4.761-08 

1.125-08 

9.153-08 

1.002-09 

7.134-11 

1.515-07 

600 

4.009-08 

9.517-09 

7.506-08 

7.725-10 

5.190-11 

1.255-07 

800 

3.549-08 

8.517-09 

6.510-08 

6.307-10 

4.050-11 

1  .098-07 

1000 

3.223-08 

7.853-09 

5.813-08 

5.319-10 

3.329-11 

9 .878-08 

2000 

2.307-08 

6.255-09 

3.937-08 

2.923-10 

2.316-11 

6.901-08 

4000 

1.470-08 

4.973-09 

2.394-08 

1.749-10 

2.877-11  1 

4.382-08 

6000 

1.051-08 

4.159-09 

1.680-08 

1.520-10 

3.249-11 

3.166-08 

8000 

8.494-09 

3.539-09 

1.288-08 

1.422-10 

3.364-11 

2.509-08 

10000 

6.964-09 

3.041-09 

1.034-08 

1.334-10 

3.354-11 

2.051-08 

TABLE  3 


CALCULATED  CROSS-SECTION  FOR  DISSOCIATIVE-RECOMBINATION,  cm2 
e  +  NO+  (v  =  0)  -*■  N  +  0 


£ 

o  2e+ 

o  2JI  III 

C  B'  2A 

2. 

0  <P 

o  Total 

.01 

1.758-14 

8.603-15 

3.092-16 

8.721-16 

2.586-17 

2.739-14 

.04 

4.599-15 

2.260-15 

8.454-17 

2.146-16 

6.934-18 

7.165-15 

.07 

3.669-15 

1.326-15 

5.498-17 

1.204-16 

4.269-18 

5.175-15 

-10 

2.537-15 

9.219-16 

4.276-17 

8.545-17 

3.170-18 

3.590-15 

.15 

1.367-15 

6.185-16 

3.414-17 

5.849-17 

2.340-18 

2.081-15 

.20 

1.302-15 

4.990-16 

2.977-17 

4.280-17 

1 

1.929-18  ; 

1.875-15 

.25 

1.032-15 

3.888-16 

2.774-17 

3.442-17 

1.710-18 

1.484-15 

.30 

9.440-16 

3.372-16 

2.690-17 

2.874-17 

1.568-18 

1.338-15 

.35 

7.302-16 

3.006-16 

2.711-17 

2.410-17 

1.465-18 

1 

1.083-15 

.40 

7.774-16 

2.558-16 

2.659-17 

2.114-17 

1.417-18 

1.082-15 

.45 

6.236-16 

2.376-16 

2.685-17 

1.860-17 

i 

1.403-18 

1 

9.080-16 

.50 

5.899-16 

2.150-16 

2.694-17 

1.639-17 

1.409-18 

8.496-16 

.60 

4.802-16 

1.810-16 

2.814-17 

1.333-17 

1.419-18 

7.042-16 

.70 

4.328-16 

1.499-16 

3.216-17 

1.086-17 

1.497-18 

6.272-16 

.80 

3.943-16 

1.322-16 

3.500-17 

8.968-18 

1.573-18 

5.721-16 

.90 

3.235-16 

1.181-16 

3.741-17 

7.279-18 

1.721-18 

4.880-16 

1.00 

2.777-16 

1.021-16 

3.785-17 

5.872-18 

1.946-18 

4.255-16 

1.20 

2.118-16 

8.498-17 

3.968-17 

3.514-18 

2.461-18 

3.425-16 

1.50 

1.357-16 

6.356-17 

5.104-17 

1.428-18 

2.005-18 

2.537-16 

2.00 

6.601-17 

3.983-17 

8.066-17 

3.073-19 

7.722-19 

1.876-16 

2.50 

3.234-17 

2.612-17 

9.338-17 

6.890-20 

1.828-19 

1.521-16 
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TABLE  4 


2  + 

DIABATIC  l  POTENTIAL  ENERGY  CURVES  FOR  NO 


R(bohrs) 

A  2£+(a.u . ) 

I  2£+(a.u.) 

H.j  j=AH(A2I+-X 

2.8 

-128.682681 

-128.676266 

.0002266 

2.81 

-128.682579 

-128.677912 

.0005507 

2.82 

-128.682566 

-128.679408 

.0009352 

2.825 

-128.682652 

-128.680042 

.0010807 

2.83 

-128.682865 

-128.680535 

.0011716 

2.84 

-128.680985 

-128.683787 

.0012275 

2.85 

-128.680987 

-128.685105 

.0012124 

2.86 

-128.680831 

-128.686532 

.0010722 

2.9 

-128.679802 

-128.692210 

.0000281 

TABLE  5 


EXCITATION  CROSS-SECTIONS  FOR  lr  STATES  OF  NO 


Channel  11  =  N  [2D]  +  0  (3P);  Channel  22  =  N  PS]  +  0  [  P] 


Scattering  Cross-Section  Matrix, 


Collision  Energy  (eV) 


o  (cra^l 


0.3375 

2.98 

1.9-07 

8.34-17 

5.3-24 

1.5-06 

5.97 

4.2-23 

1.67-16 

0.6096 

1.10 

7.6-08 

3.08-17 

2.1-24 

3.7-06 

5.97 

1.0-22 

1.67-16 

0.8817 

0.90 

5.6-07 

2.52-17 

1.6-23 

2.1-06 

5.46 

5.9-23 

1.53-16 

1.1538 

1.96 

2.9-07 

5.49-17 

8.1-24 

8.9-07 

5.24 

2.5-23 

1.47-16 

1.4259 

2.71 

4.9-06 

7.59-17 

1.4-22 

1.3-05 

6.94 

3.6-22 

1.94-16 

1.6980 

2.59 

1.4-06 

7.25-17 

3.9-23 

3.4-06 

6.31 

9.5-23 

1.77-16 

1.9702 

2.39 

1.8-07 

6.69-17 

5.0-24 

4.0-07 

6.05 

1.1-23 

1.69-16 
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TABLE  6 


MOLECULAR  CORRELATION  DIAGRAM  FOR 
0+  +  N2  •*  NO+  +  N 


o+  [4su]  +  n2 

E  =  0.0  eV 

o  [3Pg]  +  n2+  [x2xg+l 

E  =  1.94  eV 

0  [3Pg]  4-  n2+  U2nJ 
E  =  3.05  eV 

o+  [2du]  +  n2  [Xxig+] 
E  =  3.32  eV 

o  [\]+  n2+  [x2>:g+] 

E  «  3.91  eV 
o+  I2pu]  +  n2  [xxJ:g+l 

E  =  5.02  eV 

0  irlD  ]  +  N2+  lA2Eu] 

E  =  5.03  eV 

0  [3Pgl  +  n2+  [b2iu+1 

E  =  5.11  eV 

0  [3S  1  +  n2+  [X2Ig+l 
F.  =  6.13  eV 

0+  f4su]  +  N2  fA3Eu+l 
F.  =  6.16  eV 


NO+[xV"]  + 

E  =  -1 

N0+  [X1Z+1 
E  =  1. 

NO+  [  X1!"1"  1 

E  =  2. 

NO  [  X2II  1  + 
E  =  4. 

no+  fa3;:+] 
E  =  5. 

no  [x2n]  + 

F.  =  6, 

NO+  Ib3n]  - 

F.  =  6 

NO+  [w3A  1  - 

E  =  6 

NO+  [  b’  3)'+ 

E  =  7 

NO+  f  A1 X  ] 

F.  =  7 


N  [4Sul 
.09  eV 

+  N  [2Du] 
29  eV 

+  N  [ 2PU  ] 
49  eV 

N+  f3P  1 
,17  eV 

+  N  [\] 
,41  eV 

N+  [1D  ] 
g 

.07  eV 

f  N  l\] 
.14  eV 

F  N  f4S0] 
.54  eV 

]  +  N  [4SU] 
.24  eV 

+  N  [ 4Su  ] 
.44  eV 
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Atmospheric  Sciences  Laboratory 
U.  S.  Army  Electronics  Command 
White  Sands  Missile  Range,  NM  88002 

2  cys  ATTN:  DELAS-EO 

F.  E.  Niles 
M.  Heaps 


Director 

U.  S.  Army  Ballistic  Research  Laboratories 
Aberdeen  Proving  Grounds,  MD  21005 

1  cy  ATTN:  George  K.  Keller 
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Department  of  the  Army  (Cont'd) 

Commander 

U.  S.  Army  Foreign  Sciences  and  Technology  Center 
220  7th  Street,  N.  E. 

Charlottesville,  VA  22901 

1  cy  ATTN:  Robert  Jones 

U.S.  Army  Research  Office 
Durham,  NC  2770!) 


2  cys  ATTN:  R.  Mace 

Department  of  the  Navy 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  D.  C,  20350 

1  cy  ATTN:  NOP  98 5 
I  cy  ATTN:  NOP  094H 


Office  of  Naval  Research 
Department  of  the  Navy 
BOO  North  Quincy  Street 
Arlington,  VA  22217 


l  cy  ATTN:  Code  465,  K.  C .  Joiner 
1  cy  ATTN:  Code  427,  H.  Mul laney 


Commande  r 
Naval  Electronic 
Department  of  the 
Washington,  D.  C. 

I  cy  ATTN: 
l  cy  ATTN: 
l  cy  ATTN: 
l  cy  ATTN: 

1  cy  ATTN: 


Systems  Command 
Navy 
20360 

PME-117 
PME-117T 
PMK-1 17-2 1 
PMh-l 17-21A 
PMK-1 17-22 
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Department  of  the  Navy  (Cont'd) 

Director 

Naval  Ocean  Systems  Center 
Electromagnetic  Propagation  Division 
271  Catalina  Boulevard 
San  Diego,  CA  92152 

1  cy  ATTN:  Code  2200,  W.  F.  Moler 
1  cy  ATTN:  Code  2200,  Ilan  Rothmuller 
1  cy  ATTN:  Code  2200,  John  Bickel 

Naval  Postgraduate  School 
Monterey,  CA  93940 

1  cy  ATTN:  Code  0142 

1  cy  ATTN:  Code  1424 

Office  of  Naval  Research,  Dept,  of  the  Navy 
Arlington,  VA  22217 

1  cy  ATTN:  Code  421,  B.  R.  Junker 

1  cy  ATTN:  Code  465,  R.  G.  Joiner 

1  cy  ATTN:  Code  427,  H.  Mullaney 

Director 

Naval  Research  Laboratory 
4555  Overlook  Avenue,  S.  W. 

Washington,  D.  C.  20375 


1  cy  ATTN: 

1  cy  ATTN: 

2  cys  ATTN: 
1  cy  ATTN: 


Code  7700,  Timothy  P.  Coffey 
Code  7709,  Wahab  Ali 
Code  7750,  John  Davis 
Code  2627 


Commander 

Naval  Surface  Weapons  Center  (White  Oak) 
Silver  Spring,  MD  20910 

1  cy  ATTN:  Technical  Library 


Office  of  Naval  Research  Branch  Office 
1030  East  Green  Street 
Pasadena,  CA  91106 

l  cy 
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Department  of  the  Air  Force 
Commander 

Air  Force  Geophysical  Laboratory,  AFSC 
L.  G.  Hanscom  Air  Force  Base,  MA  01731 


1 

cy 

ATTN 

LKD, 

W. 

Swider 

1 

cy 

ATTN 

LKB, 

K. 

Champion 

1 

cy 

ATTN 

PHG , 

F. 

R.  Innes 

1 

cy 

ATTN 

OPR, 

R. 

A.  Armstrong 

1 

cy 

ATTN 

OPR, 

H. 

Gardiner 

1 

cy 

ATTN 

LK0, 

R. 

E.  Hoffman 

1 

cy 

ATTN 

OPR, 

A. 

T.  Stair 

1 

cy 

ATTN 

PHI, 

J. 

R.  Jasperse 

Director 

Air  Force  Technical  Applications  Center 
Patrick  Air  Force  Base,  FL  32920 

cy  ATTN:  TD 

i  cy  ATTN:  HQ  1035th  TCHOG/TFS 

Commander  AF  Weapons  Laboratory 
Kirtland  AFB,  Albuquerque,  NM  87117 

1  cy  ATTN:  H.  0.  Dogliani 

1  cy  ATTN:  J.  I.  Generosa 

1  cy  ATTN:  SUL 

Department  of  Defense  Contractors 

General  Fleet ric  Company 

TEMPO  -  Center  for  Advanced  Studies 

81b  State  Street 

Santa  Barbara,  CA  93102 

1  cy  ATTN:  T.  L.  Stephens 

1  cy  ATTN:  Warren  S.  Knapp 

1  cy  ATTN:  DASIAC 

Lockheed  Missiles  and  Space  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

1  cy  ATTN:  J.  B.  Reagan 
l  cy  ATTN:  W.  Imhof 
1  cy  ATTN:  Martin  Walt 
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Department  of  Defense  Contractors  (Cont'd) 

Mission  Research  Corporation 

735  State  Street 

Santa  Barbara,  CA  9J1U1 

1  cy  ATTN:  M.  Scheibe 
1  cy  ATTN:  D.  Sowie 

Pacific-Sierra  Research  Corporation 
1456  Cloverfield  Boulevard 
Santa  Monica,  CA  90404 

1  cy  ATTN:  E.  C.  Field 

Pennsylvania  State  University 
Ionospheric  Research  Laboratory 
College  of  Engineering 
31B  Electrical  Engineering  -  East  Wing 
University  Park,  PA  16802 

1  cy  ATTN:  John  S.  Nisbet 
1  cy  ATTN:  Les  Hale 
1  cy  ATTN:  A.  J.  Ferraro 
1  cy  ATTN:  H.  S.  Lee 

R&D  Associates 
4640  Admiralty  Way 
Marina  Del  Rey ,  CA  90291 

1  cy  ATTN:  R.  Lelevier 
1  cy  ATTN:  F.  Gilmore 
1  cy  ATTN:  R.  Turco 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  CA  90406 

1  cy  ATTN:  Cullen  Crain 

Professor  Chalmers  F.  Sechrist 
155  Electrical  Engineering  Building 
University  of  Illinois 
Urbana,  IL  61801 

1  cy  ATTN:  C.  Sechrist 
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Department  of  Defense  Contractors  (Cont'd) 

SRI  International 
133  Ravenswood  Avenue 
Menlo  Park,  CA  94023 

1  ev  ATTN:  Allen  M.  Peterson 
1  cv  ATTN:  Rav  L.  Leadabrand 
l  cv  ATTN:  J.  R.  Peterson 
l  cy  ATTN:  F.  T.  Smith 

University  of  Pittsburgh 
Dept,  of  Phvsics  and  AsLronomy 
Pittsburgh,  PA  15260 

l  cy  ATTN:  M.  A.  Biondi 
1  cv  ATTN:  F.  Kaufman 
l  cy  ATTN:  W.  Fite 
1  cy  ATTN:  .1 .  N.  Bards  ley 
l  cy  ATTN:  F..  Zipf 

General  Electric  Co.,  Space  Division,  VKSC 
Goddard  Blvd.,  King  of  Prussia 
P.  0.  Box  8555 
Philadelphia,  PA  19101 

1  cy  ATTN:  M.  11.  Bortner 

National  Oceanic  f<  Atmospheric  Admin. 
Environmental  Research  Laboratories 
Dept,  oi  Commerce 
Bon  1 der ,  CO  80302 

l  cv  ATTN:  K.  F. .  Ferguson 
1  cv  ATTN:  F.  Fehsenfeld 

Photomet r l cs ,  Inc. 

4  Arrow  Drive 
Woburn,  MA  01801 

l  cv  ATTN:  I.  L.  Kofskv 
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Department  of  Defense  Contractors  (Cont'd) 

Harvard  College  Observatory 
Cambridge,  MA  02138 

I  cy  ATTN:  A.  Dalgarno 
1  cy  ATTN:  D.  E,  Freeman 

Utah  State  University 
Logan,  UT  84322 

1  cy  ATTN:  K.  D.  Baker 

United  Technologies  Research  Center 
East  Hartford,  CT  06108 

1  cy  ATTN:  G.  A.  Peterson 
1  cy  ATTN:  H.  H.  Michels 

Stewart  Radiance  Laboratories 
1  De  Angelo  Drive 
Bedford,  MA  01730 

1  cy  ATTN:  R.  D.  Sharma 
1  cy  ATTN:  J.  C.  Ulwick 

Aerodyne  Research,  Inc. 

Crosby  Drive 
Bedford,  MA  01730 

1  cy  ATTN:  M,  Camac 

Visidyne,  Inc. 

NW  Industrial  Park 
Rurlington,  MA  01803 

1  cy  ATTN:  T.  C.  Degges 
1  cy  ATTN:  C.  H.  Humphrey 
1  cy  ATTN:  H.  J.  R.  Smith 

Physical  Sciences,  Inc. 

Lakeside  Office  Park 
Wakefield,  MA  01880 

l  cy  ATTN:  G.  Caledonia 
1  cy  ATTN:  T.  Rawlins 
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Department  of  Defense  Contractors  (Cont'd) 

Jet  Propulsion  Laboratory,  CIT 
Pasadena,  CA  91103 

1  cy  ATTN:  S.  Prasad 

BC  Space  Data  Analysis  Laboratory 
BB5  Centre  Street 
Newton,  MA  02159 


1  cy  ATTN: 

R&D  Associates 
1401  Wilson 
Ar  1  ington , 

1  cy  ATTN: 

Lockheed 
Huntsville,  AL 

1  cy  ATTN: 
1  cy  ATTN: 


S.  Guberman 

H.  Mitchell 

33B0  3 

D.  Divis 
C.  M.  Bowden 


Boulevard 
VA  22209 


Spectral  Sciences  Inc. 
99  S.  Bedford  St  reet 
Burlington.  MA,  01809 

1  cy  ATTN:  F.  Bien 
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